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ABSTRACT 

We consider the compatibility of the equations of motion which follow from d'Alembert's 
principle in the case of a general autonomous non-holonomic mechanical system in dimen- 
sions, with those equations which follow for the same system by assuming the validity of a 
specific variational action principle, in which the non-holonomic conditions are implemented by 
means of the multiplication rule in the calculus of variations. The equations of motion which 
follow from the principle of d'Alembert are not identical to the equations which follow from 
the variational action principle. We give a proof that the solutions to the equations of motion 
which follow from dAlembert's principle do not in general satisfy the equations which follow 
from the action principle with non-holonomic constraints. Thus the principle of d'Alembert and 
the minimal action principle involving the multiplication rule are not compatible in the case 
of systems with non-holonomic constraints. For simplicity the proof is given for autonomous 
systems only, with one general non-holonomic constraint, which is linear in the generalized 
velocities of the system. 
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1 Introduction 



Hamilton's principle for mechanical systems with non-holonomic constraints has recently been 
discussed by Flannery [1]. In particular a variational formulation of the equations of motion 
of a mechanical system was discussed both for holonomic and non-holonomic constraints. It 
was shown that while the equations of motion for a system with holonomic constraints can 
be obtained as variational equations, with the constraints being taken into account by the 
multiplication rule in the calculus of variations [2], the corresponding procedure with non- 
holonomic constraints leads to equations which differ from the equations of motion which follow 
from the well-known principle of d'Alembert. 

The problems discussed by Flannery are not new; they have been discussed in the literature 
at least since Hertz's textbook [3], in which the use of variational principles in mechanics was 
questioned. We refer in particular to an early paper by Holder [1], in which essential differences 
between holonomic and non-holonomic systems was discussed. Two later papers published by 
Jeffreys [5] and Pars [0] considered again Hamilton's principle for non-holonomic systems, and 
proposed rectification of previous papers in which the variational procedure (action principle) 
involving the multiplication rule had been proposed for systems with non-holonomic constraints. 

The papers by Flannery, Jeffreys and Pars referred to above, contain several references to 
papers in which one has advocated the use of a variational action principle involving the mul- 
tiplication rule for non-holonomic systems. We should also mention a paper by Berezin [7|, in 
which it is taken for granted that the same minimal action principle is valid both for holonomic 
and non-holonomic systems, with the constraints taken into account using the multiplication 
rule in the calculus of variations. 

Even though the equations of motion following from the principle of d'Alembert and from 
the variational action principle with non-holonomic constraints are different in form, one may 
still argue that the equations in question may have the same solutions. It was demonstrated 
by Pars [S], that the solutions to the equations in question differ from one another, in the 
case of a simple example of a non-holonomic system in three dimensional space. As such, this 
validates the assertion that the solutions to the two different types of equations of motion are 
different in general, at least in a space of three dimensions. To the best of our knowledge it 
has not been proved for general systems that the the two kinds of equations of motion referred 
to above in general have different solutions. We will supply a proof of this fact in this paper, 
for a general autonomous system with a finite number of degrees of freedom, restricted only 
by reasonable smoothness conditions. For simplicity we consider in our proof only the case of 
one non-holonomic constraint, which is taken to be linear and homogeneous in the generalized 
velocities of the system. 

The problem we consider ought to be of interest in Classical Mechanics, but may also be 
of interest for the question of quantization of non-holonomic systems. One standard method 
of quantization of a classical system, namely the functional integral method, relies explicitly 
on the existence of an action principle for the classical system in question. This is also true 
for the case of canonical quantization, which requires the existence of canonical coordinates 
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and momenta and a proper Hamiltonian. This, in turn, requirers the existence of a proper 
Lagrangian, and hence the existence of an action principle. 

The proof given in this paper imphes definitely that the standard action principle, which is 
valid for unconstrained systems, and which can be generalized to cover the case of holonomic 
constraints by using the multiplication rule in the calculus of variations, is not consistent with 
the principle of d'Alembert if the non-holonomic constraints are implemented by means of 
the multiplication rule. Hence this action principle can not be taken as a starting point for 
quantization of non-holonomic systems, unless one dismisses the principle of d'Alembert in the 
case of non-holonomic systems, and assumes the validity of the action principle involving the 
multiplication rule instead. 

In the next section we consider the Lagrange equations of motion which follow from the 
principle of d'Alembert for an autonomous mechanical system with both holonomic and non- 
holonomic constraints. We also derive the corresponding variational equations, which are ob- 
tained by using the multiplication rule for the implementation of the constraints. 



2 Lagrange equations with constraints 

Consider an autonomous mechanical system with independent generalized coordinates q = 
{q^, g^), and velocities q = {q^, g^). We denote the kinetic energy of the system by T, and 
the generalized applied forces on the system by Qa, A = 1, N. The principle of d'Alembert 
(see e.g. the classical texts by Goldstein [8j or Whittaker [9j) then gives the following equation. 



N 

-f^\dt \ dq^ J dq^ 



where the quantities Sq^ are virtual displacements of the system. If the virtual displacements 
5q^,A = 1,...,N, are independent, then the Eqn. ([1]) results in the Lagrange equations of 
motion. 

We generalize to systems with 1 < m < independent non-holonomic constraints, which 
are taken to be linear and homogeneous in the velocities. The non-holonomic constraints are 
thus of the following form, 

N 

E«a(?)^^ = 0, ^ = l,...,m<iV, (3) 

A=l 

where the quantities a\{q) are given functions of the variables q = {q^, q^). 

The derivation given below of the equations of motion for this non-holonomic system can 
be found in the textbook by Whittaker [9]. 
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Implement the constraints ([3]) by regarding the system to be acted on by both the external 
applied forces Qa and by certain additional forces of constraint Q^, A= 1, . . . , N, which force 
the system to satisfy the non-holonomic conditions ([3]). Equation ([1]) is then replaced by the 
following equation, 

" d f dT\ _ dT 

In Eqn. (jl]) the virtual displacements Sq"^, A = 1, . . . ,N, can now be regarded as independent. 
Thus one obtains the equations of motion, 



E^ItIjSS -JC3-<3^-% V = 0, (4) 



The forces of constraint, Q'^, A = 1, . . . ,N, are a priori unknown, but they are such that, in 
any instantaneous displacement Sq'^jA = 1,...,N, consistent with the constraints ([3]), they 
do no work. The non-holonomic constraints ([3]) imply the following conditions on the possible 
instantaneous displacements 6q^, A = 1, . . . , N, 

N 

E«a(9)V = 0, ^ = l,...,m<iV. (6) 

A=l 

For any instantaneous displacements Sq'^yA = 1, . . . ,N, which satisfy the conditions ([6]), the 
work 6W' done by the constraint forces Q^, A = 1, . . . , is equal to zero, i.e. 

5W' ^ £ Q'^5q^ = 0. (7) 

A=l 

The conditions ([6]) and ([7]) together imply that 

m 

Q'A = T.^^^Aiq),A = l,...,N, (8) 

i=l 

where the quantities Aj, z = 1, . . . , m, are time- dependent parameters. Eqns. ([5]) have thus been 
reduced to the following form. 

These equations of motion are consequences of the principle of d'Alembert. One should 
still add the m equations of constraint ([3]) to the equations of motion above. There are thus 
altogether N + m equations for the determination of iV + m quantities q^(t),A = 1, . . . , N, and 
Aj(t), i = 1, . . . , m, when appropriate initial conditions for the quantities g^, and q^, ...,q^ 
are given. 

One should observe that in the argument above, that it is not required that the constraint 
equations ([3]) be in force under general variations qj — > qj + Sqj; the constraints ([3]) are only 
imposed on the actual motion of the system. 
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It is convenient to assume that the external apphed forces can be derived from a generahzed 
potential This entails unimportant loss of generality. In what follows we thus assume the 
existence of a potential V such that 

Using the notation 

Loiq,q)^T-V, (11) 
we then rewrite the equations as follows, 

d ( dLo{q,q)\ dLo{q,q) 



dt \ dq^ I dq^ 



^A,a^((/), A=1,...,N. (12) 



i=l 



The quantity LQ{q,q) defined above in Eqn. flTTl) depends on the kinetic energy and on the 
external applied forces of the system under consideration. If there were no non-holonomic 
constraints then the quantity Lo{q, q) would be the Lagrangian of the system. 

It should be noted, that from he equations (fT^ and the constraints follows a first integral 
for the system. 



N 



.A dLo{q,q) 



where E is the constant energy of the system. 

We now consider the variational action principle discussed above, in which the non-holonomic 
conditions ([3]) are taken into account by means of the multiplication rule in the calculus of vari- 
ations. Consider the following action functional, 

5o:= JdtLo{q,q). (14) 

The action principle referred to previously is simply the requirement that the action functional 
( fT4l) be stationary when the non-holonomic constraints ([3]) are in force throughout a suitable 
region D in configuration space, i. e. when the conditions 

N 

E «A(9)g^ = 0, ^ = 1, m < TV ; (g\ g^) G D, (15) 
A=l 

are in force. Using the multiplication rule in the calculus of variations, the action principle 
formulated above becomes equivalent to the following free variational problem. 



m N 



i=l A=l 



6 J dt 

which involves m Lagrange multipliers iJ,i,i = 1 , . . . , m 



0, (16) 
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The variational equations following from Eqn. f|T6|) are, 



d / dLo{q,q) \ dLo{q,q) _ 

dt [ dq^ J dq^ ^ ^ 



B=l 



QqB Qq/ 



The N equations f|T7|) together with the m conditions f|T5|) are supposed to determine the 
quantities q^, ...,q^ and the Lagrange multiphers /xi, . . . , fim, when appropriate initial conditions 
for q^, ...,q^ and g^, are given. 



The variational equations f[T7j) are not identical to the d'Alembertian equations of motion 
([T2|) . The assumptions underlying the variational equations and the d'Alembertian equations 
of motion, respectively, are also basically different. The non-holonomic constraints are only 
supposed to be valid for the actual motion when one applies the principle of d'Alembert to 
derive the equations of motion, whereas these non-holonomic conditions are supposed to be 
valid throughout a whole appropriate region D when one considers the variational problem 
(|T6|) . as indicated in Eqn. (fTSll above. 



It is readily seen that if the system under consideration is holonomic, i.e. if the one-forms 
occurring in Eqn. ([6]) are integrable, then the d'Alembertian equations (fT2|) and the variational 
equations (fT7|) are equivalent. One reaches the same conclusion when the one-forms in ([6]) can 
be made integrable by multiplication with integrating factors. We will return to the question 
of existence of such integrating factors in the next section. The forms are integrable if the 
following conditions are satisfied. 

When the conditions ( fTSl) are in force, then the d'Alembertian equations of motion ( fT2l) and the 
variational equations ( ITTI) become identical upon a change of notation: /ij — Aj, i = 1, . . . , m. 

The variational problem f|T6l) is readily recognized to be a variational problem with holo- 
nomic constraints when the integrability conditions ( fT8|) are in force, since then there exist 
functions \E'*(g) such that 

d_ 

dq^ 

Here and in what follows we tacitly assume that the configuration space is a contractible man- 
ifold. The non-holonomic constraints ([3]) then become equivalent to the following conditions. 



<iQ) = 7r^'^\Q), ^ = l,...,m. (19) 



j^^\q) = 0, 2 = l,...,m, (20) 

from which the following holonomic constraints follow, 

^\q) = C\ z = l,...,m, (21) 
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where the quantities are constants. From the equations f|T9l) and fl2T|) follows that the 
variational problem f|T6|) is equivalent to the following problem, with holonomic constraints 
implemented by means of the multiplication rule, 



6 / dt 



0. 



(22) 



The use of the principle of d'Alembert and the variational action principle with non- 
holonomic constraints give rise to basically different problems, which also result in different 
equations of motion. When the constraints are holonomic, both problems yield identical equa- 
tions of motion, as was demonstrated above. 

It is legitimate to ask whether the d'Alembertian equations of motion and the variational 
equations, respectively, in the non-hoholonomic case nevertheless can have the same solutions 
in general. It was shown by Pars [B] in the case of a specific simple example in three- 
dimensional configuration space, involving one non- holonomic constraint, that the solutions 
of the d'Alembertian equations of motion and the corresponding variational equations do not 
coincide in general. In the next section we consider the case of a general A^-dimensional system 
satisfying the conditions outlined in the Introduction, with one non- holonomic constraint, and 
prove that the general solutions to the d'Alembertian and variational equations, respectively, 
do not coincide for any > 3. 



3 Nonequivalence of the principle of d'Alembert and the 
variational action principle 

We now consider the d'Alembertian equations of motion (fT2l) and the variational equations 
(fT7|) . respectively, but with only one non-holonomic constraint, 

N 

J2 aAiq)q^ = 0, (23) 

A=l 

where the functions aA{q), A = 1, . . . , N, are arbitrary functions of the variables q = [q^, . . . , g^), 
except for the condition that not all of the quantities a^(g) vanish identically. Naturally it is 
also assumed the the functions aA{q),A= 1, . . . ,N satisfy appropriate smoothness conditions. 
The equations (fT2l) then reduce to the following set of d'Alembertian equations of motion, 

d (dLo{q,q)\ dLo{q,q) \ . . .r foA\ 

where A is a time dependent parameter. 

It is convenient to introduce the following notation, 

, . , ^ dttAjq) daBjq) . „ . 
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The variational equations ([T7|l reduce to the following set of variational equations in the case 
of one non-holonomic constraint, 



d ( dLo{q,q) \ dLo{q,q) . ^ 

^ j ^^^ = /ia^(g)+/i|_^MAB(g)g = 0, A = 1, . . . , N., (26) 

where ia a Lagrange multiplier. 

The question is now whether the equations fl24l) and (126|) . respectively, can have the same 
solutions for . . . , in general, despite the fact that these equations are not identical. 
It will be shown below that this is not the case. 

Naturally, the functions L(){q, q) and ayi(g) entering into the equations and (126|) . respec- 
tively, will have to satisfy appropriate smoothness conditions in order that these equations may 
have solutions. We will not enter into a discussion of such smoothness conditions, but rather as- 
sume that the equations fl2^ and fl26|) . respectively, have e.g. C^-solutions in some appropriate 
time-interval, for given initial values for the coordinates q^ and the velocities q^,A= 1, . . . ,N, 

= ^0, = go^, ^ = 1, . . . , iv. (27) 



t=to 



t=to 



The initial values q^ and q^ at t = to are free parameters within an appropriate region of the 
configuration- and velocity space, except for the restriction 

N 

E (^A{qo)q^ = 0. (28) 

A=l 

The condition (128|) is a consequence of the non-holonomic constraint (123|) . 

In what follows we use the phrase "general initial value conditions" for the conditions which 
have been described and defined above in Eqns. ( 1271) and ( l28l) . 

If the functions aA{q) in Eqn. fl23l) satisfy the conditions 

MAB{q)=0, A,B = 1,...,N, (29) 

for all q in an appropriate contractible domain D, then the condition fl23l) is integrable, i.e. 
there exists a function \E'(g) ^0 such that 

= E aA{q)q^ = |vl/(g(t)) ^ vl/(g) = C, (30) 

where C is a constant. In this case the system is in fact holonomic. If the integrability condition 
(!29|) is not in force, there may exist an integrating factor ^{q) ^0, such that the condition 

N 

E $(g)aA(g)g^ = 0, (31) 

A=l 
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is integrable. In this case the system is reducible to a system with a holonomic constraint 
simply by replacing the constraint fl23|) by the equivalent constraint fl3T|) . i.e. by making the 
replacement 

aAiq) ^ Hq)aAiq) (32) 

throughout. 

We notice in passing that the non-holonomic constraint (l23l) in a two-dimensional configura- 
tion space is always reducible to an equivalent holonomic constraint by means of an integrating 
factor. In what follows we thus only consider the cases = 3, 4, . . .. 

The necessary and sufficient conditions for the existence of an integrating factor ^{q) are 
the following conditions [lOj, 

aA{q) MBciq) + asiq) McA{q) + ac{q) MAsiq) = 0, A, B,C = 1, . . . , N. (33) 

The integrability conditions formulated above are neatly expressed using exterior derivatives 
and differential forms, We refer to e.g. the text-books by Flanders [Ll\ or Fleming [12] for an 
exposition of the appropriate formalism. 

After these preliminaries we will now prove the following: 

Theorem of non- equivalence: If the constraints (123!) are truly non-holonomic, i.e. neither in- 
tegrable nor reducible to integrable constraints by means of an integrating factor, then the 
d'Alembertian equations (1241) and variational equations fl26p do not have coinciding solutions 
for arbitrary general initial value conditions fl27j) and fl28|) . The dimensionality A^ of configura- 
tion space is supposed to be a finite integer, satisfying the condition A^ > 3. 

We assume that the d'Alembertian equations have unique smooth {e.g. C^) solutions 
{q^{t), . . . ,q^{t)) in an appropriate time-interval, satisfying the non-holonomic constraint fl23|) 
and the general initial value conditions fl27|) and fl28l) . We call such solutions general solutions 
in what follows. 

The method of proof is by reductio ad absurdum, i.e. we make the assumption that the 
variational equations fl26|l have solutions which coincide with the general solutions of the 
d'Alembertian equations flMl) . and show that this assumption leads to contradictions. 

Assume that the equations fl26l) and fl2^ have coincident solutions. By subtracting Eqn. 
fl26|) from Eqn. fl24|) . one obtains the following equations, 

N 

{X-ii)aA = li J2MABiq)q'' = 0, A=1,...,N, (34) 

B=l 

which have to be satisfied by the general solutions {q^, . . . , q^) of the d'Alembertian equations 

Let us first note that one must necessarily have /z ^ in the equations (IMI) above, since 
otherwise one would have 

XaA = 0, A = l....,N. (35) 
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The condition (l35l) above implies that A = 0, which can not be true in the case of a general 
non-holonomic constraint. Hence /i ^ 0. 

We introduce the notation 

r ;= (36) 
The conditions (IMl) are then equivalent to the following equations, 

N 

J2MABiq)q'' = raA, A=1,...,N, (37) 

B=l 

which is a set of linear algebraic equations in the variables q^,A= 1, . . . , N. 

The solutions q^,A = 1, . . . , N, of the algebraic equations (1571) can be obtained in a fairly 
explicit form, by using standard properties ^13j of skew (antisymmetric) matrices. Below we 
give a brief account of the properties needed in our analysis of the equations (l37|l . A more 
detailed account is given in Appendix A at the end of this paper. 

Let M ^ be an arbitrary skew N x N matrix, with real valued matrix elements Mab, 
labeled by integers A, B in the range (1, . . . , A^), where iV > 3. Since M is skew, its rank is 
necessarily an even integer 2p < N, Since M ^ we must also have p > 1. Consider the matrix 

N 

{M^)ab ■■= E MacMcb. (38) 
c=i 

The matrix is real and symmetric, and hence it has eigenvalues —nl, with corresponding 
ortho- normal eigenvectors e^, v = 1, . . . , N. The eigenvalues —k"^ are non-positive, since 
is also negative semidefinite. The numbers Hi, are therefore real, and can be chosen to be 
non-negative without loss of generality. The non- negative quantities Kj,,z/ = 1,...,A^, are 
determined by the roots of the equation 

det{M^ + k11) = 0. (39) 

The matrix has the same rank as the matrix M . Hence, 

> 0, // = l,...,2p, (40) 

and, if 2p < iV, 

K^ = 0, z/ = 2p+l,...,A^. (41) 

As shown in Appendix A, one can form a set of ortho-normal basis vectors h^^v = 1, . . . ,N, 
which satisfy the following equations, 

N N 

E MABb2u-l,B = HuhuA'^ E ^ABhu^B = -I^ub2u-1,A, 1^ =l,...,P, (42) 
B=l B=l 



and, if 2p < N, 



N 



Y,MABbu,B = 0, u = 2p+l,...,N. (43) 



B=l 
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The basis vectors are ortho- normalized in the following inner product, 

N 



{b„ K) := E = := I J (44) 



A=l 



The basis- vectors b^,iy = 1, . . . ,N are determined by the elements Mab of the matrix M. Thus 
the components = 1, . . . , N; A = 1, . . . ,N, are in principle known quantities once the 

matrix elements Mab, A,B= 1, . . . , N are given. 

We now return to the equation fl57|) . which was shown above to be a necessary consequence 
of the assumption that Eqns. fl24|) and fl26|) . respectively, have coincident general solutions. 

It is convenient to consider separately the case F = and F ^ 0. 

The case F = 0. 

When F = 0, Eqns. (1571) read as follows, 

N 

J2MAB{q)q'' = 0, A=1,...,N, (45) 

B=l 

We first consider the case when the matrix M is regular, i.e. det(M^5(g)) ^ 0. This can 
happen only if the rank 2p of M equals A^, in which case N is an even integer. The only 
solutions of Eqns. (H5|) are then the following, 

= 0, A=1,...,N. (46) 

However, the solutions P6l) are not possible, since they are not consistent with the general 
initial value conditions f l27l) and ( !28l) . 



Let then the rank 2p of M be less than N. The general solution of Eqns. fl37l) is in this case 
a linear combination of the N — 2p basis vectors b,y, u = 2p+l,...,N, i.e., 



q^= J2 ^-Ka, A=1,...,N, (47) 

u=2p+l 

where the quantities 7,^, = 2p + 1, . . . , are free parameters. 

From the orthogonality conditions flH|) follows then that the solutions 0471) are orthogonal 
to the basis vectors when /i = 1, . . . ,2p, i.e., 

N 

{b^,q)=J2b,,Aq'' = 0, t, = l,...,2p. (48) 

A=l 

The conditions fHHl) must in particular also be satisfied for the initial values of q^ and q^,A = 
1, . . . , N, at t = to, i.e.. 



A=l 



S K'^^1t=to=°' /^ = l'---'2P- (49) 
10 



Since the number 2p is in the range (2, . . . , — 1), the number of conditions in fHOj) is at most 
— 1, but at least 2. Thus the solutions fHTl) are not possible, since they involve at least two 
special conditions of the form fH9|) on the 2N initial values and Qq^A = 1, . . . , N. This is 
not consistent with the general initial value conditions (1271) and (1251) . 

We have now shown that the case F = leads to contradictions in any A^-dimensional 
configuration space of dimension > 3. It remains to consider the case F ^ 0. 

The case F ^ 0. 

We return to the equations (157|) . where F ^ 0, 

N 

Y,MAB{q)f = TaA] F ^ 0, A = 1, . . . , iV. (50) 

B=l 

We multiply the equations flSUl) above from the left with the basis vector components K^a^ 
where v < 2p, and sum over A in the range (1, . . . , A^). Using the antisymmetry of the matrix 
elements Mae and the equations (H2ll . one obtains, 

Kuib2u~i,q) = 'r{b2u,a), u = l,...,p. (51) 

and 

Ku{b2u,q) = -'r{b2u-i,a), z/ = l,...,p. (52) 
It is convenient to rewrite the conditions (jSTj) and fl52|) above as follows, 

N 

{b„ q) = J2 K,A = F e^, = 1, . . . , 2p, (53) 

A=l 

where 

p 

■= i^u^[{b2u,a)Sf,,2u-i - (b2u-i,a)6^^2u], /i=l,...,2p. (54) 

It should be observed that the expressions for the quantities 6^ defined above in Eqn. 
contain only quantities which in principle are known when the components a^, ^ = 1, • • • , A^, 
occurring in the non-holonomic constraints fl23|) are given, namely the non- vanishing quantities 
as well as the components of the basis vectors b2u-i and 621/, for u = 1, . . . ,p. Consider now 
the conditions obtained from the equations fl55]) for any two distinct values of ^i, e.g. = a 
and /J, = (3. Eliminating the unknown but non-vanishing quantity F from these equations one 
obtains the following p{2p — 1) conditions, 

(6,, q)efs - {bfs, g)e« = 0, a, /3 = 1, . . . , 2p, (55) 

which hold for all t, and hence in particular for the initial time to- 

Since p > 1, there is at least one condition of the form flSSl) involving the initial values for 
q^ and q"^, A = 1, . . . , N, ai t = Iq. Hence the solutions q"^, A = 1, . . . ,N of the equations flH?]) 
are not possible, since they involve at least one special condition of the form (!55|) on the 2A^ 
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initial values Qq and qQ,A= 1, . . . ,N, and this is not consistent with the general initial value 
conditions fl271l and fl2H|) . 

We have now shown that the also the case F ^ leads to contradictions in any A^- 
dimensional configuration space of dimension > 3. 

Conclusion of the proof. 

We have shown that the equations (1371) . which are an inescapable consequence of the as- 
sumption that the d'Alembertian equations of motion (p^ . and the variational equations (l26l) . 
respectively, have coincident general solutions, are in contradiction with the requirement that 
the initial values and at t = to for general solutions must be considered as free param- 
eters in an appropriate region, except for the restriction implied by the given non-holonomic 
constraint fl23|l at the initial time to- We have arrived at this contradiction when the constraint 
fl23|) is truly non-holonomic. The contradiction proves the assertion that the d'Alembertian 
equations and variational equations fl2B]) do not have coinciding general solutions in any 
A^-dimensional configuration space with A^ > 3, if the constraints fl25]) are truly non-holonomic, 
i.e. neither integrable nor reducible to integrable constraints by means of an integrating factor. 
This concludes the proof. 

It was demonstrated in Sec. 2 above, that the d'Alembertian equations of motion fll2l) and 
the variational equations (1T71) become identical after a simple change of notation, if the con- 
straints ([3]) are holonomic, or reducible to a set of holonomic constraints. When the equations 
in question are identical, they naturally have identical solutions, as long as the solutions exist. 
This conclusions holds a fortiori in the case of one constraint. So, we have in fact proved that 
the d'Alembertian equations of motion fl24p and the variational equations (12^ . respectively, 
have coincident general solutions in the case of one constraint of the form fl23|) if and only if 
the constraint in question is either holonomic or reducible to a holonomic constraint by means 
of an integrating factor. 

We finally make the observation that it was sufficient for the proof of the theorem above, 
to confront the consequences of the equations fH5|) and fl50|) . respectively, with the requirement 
that the general solution {q^{t), . . . ,q^{t)) of the d'Alembertian equations of motion fl24l) must 
be specified by the general initial value conditions fl27|) and (1251) . which involve 2A^ — 1 initial 
values as free parameters in an appropriate region. This requirement by no means exhausts the 
consequences of the equations (jl51) or (15U]) . respectively. 

4 Summary and discussion 

It has been proved in this paper that the d'Alembertian equations of motion for a fairly general 
system with non-holonomic constraints, in three or more dimensions, are not equivalent to the 
corresponding variational equations with the constraints implemented by the multiplication rule 
in the calculus of variations. The non-equivalence means two things, namely that the equations 
in question are not identical in form, and also that they do not have coincident general solutions. 
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The non-coincidence of the general solutions of the d'Alembertian and variational equations, 
respectively, has been demonstrated in a special case in three dimensions in a paper by Pars 
[6], but has not until now, to the best of our knowledge, been demonstrated for general systems 
in dimensions with > 3. 

The action principle, which has been proposed several times in the literature for systems with 
non-holonomic constraints, in analogy with the valid action principle for systems with holonomic 
constraints, is thus not consistent with the principle of d'Alembert. It is possible to resolve 
this problem of inconsistency by dismissing the principle of d'Alembert as a valid principle for 
systems with non-holonomic constraints. This possibility has been explicitly rejected for good 
reasons in the papers by e.g. by Jeffreys [5] and Pars [6J, which were referred to previously. 
Also e.g. in the classic text by Whittaker [9], the principle of d'Alembert is considered to be 
fundamental. 

It is possible that the d'Alembertian equations of motion in specific cases may have solutions 
which coincide with the solutions of the corresponding variational equations for some special 
initial values. In such cases the solutions in question do naturally also follow from the action 
principle involving the multiplication rule in the calculus of variation. However this action 
principle would in such a case not be universal, but would be tied to particular cases with 
special initial values. 

It is legitimate to ask whether there may nevertheless be some action principle governing 
systems with non-holonomic constraints. The analysis of this question can only use the equa- 
tions of motion in the non-holonomic case, and will then have to rely on such results which 
are available for the so-called inverse problem in the calculus of variations [H], [15]. Thus any 
variational principle which might be found in this manner for some specific system would not 
be universal. 

In a previous paper [16] we have considered the existence of a variational principle for the 
simple three-dimensional example with one non-holonomic constraint considered by Pars [B], 
and have shown that there is indeed an action principle in that case, in the sense that there 
is a related system in two dimensions which can be formulated in terms of an action principle 
involving a Lagrangian and also a Hamiltonian. This related system yields the equations of 
motion for the three-dimensional example. The Lagrangian for this case is non-unique in a 
non-trivial way, however. 
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Appendix A. Properties of skew matrices 



In this Appendix we give a brief account of those properties of skew (antisymmetric) matrices 
that are needed in Section 3 of this paper. For a detailed account of the material presented 
here, we refer to a text-book by Greub |13j . 

Let M be an arbitrary skew N x N matrix, with real- valued matrix elements Mab, labeled 
by integers A, B in the range (1, . . . , A^), where > 3. Consider the matrix M^, 

N 

{M^)ab := E MacMcb. (A.l) 
c=i 

From the anti-symmetry of the matrix M follows that the matrix is symmetric and negative 
semi-definite. Hence there exists a set of ortho-normal eigenvectors {e,^} , N = 1, . . . , A, which 
solve the following eigenvalue equations, 

N 

E iM^)ABeu,B = -f^l e.,A, z/ = 1, AT ; A = 1, 2, . . . , A, (A.2) 

B=l 

where the quantities K,y are real numbers, which without loss of generality can be assumed to 
be non-negative. The eigenvalues and eigenvectors are labeled by Greek letters such as /x, z^, . . ., 
and the single indices A, B, . . . label the components of the eigenvectors, or, more generally, 
the components of any vectors in the domain of M. The inner product {u, v) of any two such 
vectors u and v is taken to be the following Euclidean inner product, 

N 

{u,v) := E uaVa- (A.3) 

A=l 

Since one uses an Euclidean inner product it is not necessary to distinguish between upper 
and lower indices of the components of vectors in the domain of M. We will use the inner 
product (1A.3P also if one or both of the vector components are labeled by upper indices. The 
eigenvectors 6^,1^ = 1, . . . ,N, are ortho-normalized, 

("-'"^) = ^-^ = \ otherwise. ^^'^^ 

The rank of a skew matrix is always an even non-negative integer. The rank of is the 
same as the rank of M. Let this rank be denoted by 2p in what follows. Thus there are 2p 
eigenvalues — /t^ in the equations flA.2p which are different from zero, i.e. 

> 0, V = 1, ...,2p. (A.5) 

If 2j9 < A^ there are N — 2p remaining eigenvalues which are zero, i.e., 

K^ = 0, u = 2p+l, ...,N. (A.6) 
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Define a set of basis vectors b,y as follows, 

b2,i-i,A ■= e^^A, /i = 1, . . . ,p, 

TV 



A 



B=l 

l,...,iV, 



and, wlien 2p < N, 

bu,A ■= e^^A, ly = 2p+ 1, . . . ,N; A = l, 
The basis- vectors 6^, /i = 1, . . . , are ortho- normal, 

(6^,6^) = S^i,. 

We finally note that the equations (]A.7p and (lA.SP imply the following conditions. 



N 



N 



MABb2v-l,B = Hub2iy,A; MABb2u,B = ~Hub2u-l,A^ ^ = 1) ■ • • 



B=l 



B=l 



and, if 2p < N, 



N 



E MABbu,B = 0, iy = 2p+l,...,N. 



B=l 



In the basis {bf^}\ , the matrix M has the following so-called normal form. 



M 



1 ° 






. ...\ 


-Ki 


. 










Kp 








. —Kp 








. 




V ... 






. / 



which is an immediate consequence of the equations flA.lOp and flA.lip above 



(A.7) 

(A.8) 
(A.9) 

(A.IO) 
(A.ll) 



(A.12) 



Any vector v in the domain of M can be expanded in the basis {b^}i . It is convenient to 
write such an expansion as follows, 

p N_ 

- (A.13) 



^=1 ^=2p+l 



with the understanding that the sum Zl/i=2p+i(- • •) is absent if 2p = N. The quantities a^, 
and 7^ are readily obtained with the aid of the ortho-normality conditions (lA.9p . 

= {b2i.-i,v), = {b2u,v)] u = l,...,p, (A. 14) 

and, when 2p < N, 

-f, = {b,,vy, u = 2p + l,...,N. (A.15) 



We will make frequent use of all the equations (1A.1QI - lA.llI) and (]A.13I - IA.15I) above in 
the main body of the text. 
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